Edge exponent in the dynamic spin structure factor of the Yang-Gaudin model 
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The dynamic spin structure factor S(k,uj) of a system of spin-1/2 bosons is investigated at ar- 
bitrary strength of interparticle repulsion. As a function of ui it is shown to exhibit a power-law 
singularity at the threshold frequency defined by the energy of a magnon at given k. The power-law 
exponent is found exactly using a combination of the Bethe Ansatz solution and an effective field 
theory approach. 



The remarkable progress achieved by the theory of 
one-dimensional (ID) quantum fluids is rooted in the 
fact that dimensionality imposes severe constraints on 
the fluid's low energy excitation spectrum. Due to these 
constraints the investigation of the low-energy dynamics 
of the fluid reduces to choosing the effective field the- 
ory from a limited number of universality classes. Per- 
haps the most ubiquitous (and most thoroughly investi- 
gated) is the universality class called the Luttingcr Liq- 
uid Other non- trivial examples include states with 
non-abelian currents, spin-incoherent [E IE 0] and ferro- 
magnetic liquids [1, IE LE IE Q • For au sucn cases there 
exist well developed analytical methods allowing one to 
calculate infrared asymptotics of dynamical correlation 
function, spectral properties, and scaling dimensions of 
local observables. 

In a series of recent pap ers [E IE 0, IE IE EE U, 12, 
EE EE EE EE EE Hi it has been found that the 
dimensionality constraints and the resulting universality 
may extend far beyond the low energy sector of the ex- 
citation spectrum of the fluid. It was shown that there 
exist a curve ui-(k) in the (k,u>) space at which spectral 
functions exhibit power law singularities of the type 

S{k,u>) ~ c(k)9[uj - oj-(k)]{oj - uj-(k)] A{k) . (1) 

Here 9{x) is the Heaviside step function, A(k) and c{k) 
are some momentum-dependent functions, and to-(k) is 
the energy of the lowest excited state of the fluid at a 
given momentum k. In a generic ID fluid > for 

all k except a discrete set of points defined by the Lut- 
tinger theorem. The spectrum of momentum-dependent 
anomalous exponents A(fc) in Eq. ([1]) is a natural gen- 
eralization of the spectrum of scaling dimensions of the 
low-energy effective theory. Understanding the structure 
of the spectrum of A(fc) will greatly advance the theory of 
ID quantum fluids. There are several approaches to this 
problem. In Refs. [HI, EE] perturbation theory was used 
to get A(k) in a fermionic system. In Refs. IE 13. 14\ 
an effective field theory approach establishin g a link with 
the mobile quantum impurity problem pol. 1 2 ll. |22| was 
proposed. This approach was complemented by Bethe 
Ansatz (BA) calculations for several integrable models: 
Calogero-Sutherland [H|, Heisenberg [IEEE El, and 



Licb-Liniger [l8| . Constraints on A(k) implied by sym- 
metries of microscopic Hamiltonian were discussed in 
Refs. @, Ell- 
in a recent work [BJ on the dynamical properties of 
a strongly repulsive ferromagnetic Bose gas observable 
phenomena such as spin trapping and gaussian damping 
of spin waves were predicted and a link between these 
phenomena and the singular behavior, Eq. ([T]), of the 
dynamic spin structure factor was established. It was 
shown that at infinite point-like repulsion and for k — > 0, 



A(fc) ~ -1 



K ( k 



(2) 



where K is the Luttinger parameter and kp — 7rp$ with 
Po being average particle density. Assuming the validity 
of Eq. ([2]) at a large but finite repulsion, a crossover be- 
tween trapped and open regimes of spin propagation was 
characterized completely. A different approach to the dy- 
namics of the same system proposed in Ref. confirmed 
Eq. @ • The approach of Ref. Q was further developed 
in Ref. 0], demonstrating that for infinite point-like re- 
pulsion A(fc) has the form @ for arbitrary k. In Ref. Q 
the small k expansion of A(fc) was shown to have the 
form ^) for arbitrary interparticle repulsion. However, 
the case of arbitrary k and interparticle repulsion remains 
unexplored. 

In this paper we investigate the behavior of A(fc) and 
U-(k) for the dynamic spin structure factor of a ferro- 
magnetic system of spin-1/2 bosons interacting through a 
point-like repulsive potential of arbitrary strength. This 
system is described by the Yang-Gaudin model [23| . 
Combining the Bethe Ansatz with an effective field the- 
ory we obtain our main result: explicit expressions (I20[) 
(f2"2")l for A(k) at arbitrary k and interparticle repulsion. 

The Hamiltonian of the Yang-Gaudin model is 



H = 



d.r 



(3) 



where ipm)(%), ^m)( x ) are canonical Bose fields satisfy- 
ing periodic boundary conditions on a ring of circumfer- 
ence L, and p(x) is the total particle density operator. 
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We consider the dynamic spin structure factor 
S(k,uj)= [ dxdte l{uJt - kx) (1[ |s+(x,t)s_(0,0)| ft). (4) 



Here s + (x) = ip]-(x)ip[(x) is the local spin raising oper- 
ator, and S-(x) = [s + (x)]ft The average in Eq. (j4|) is 
taken with respect to a fully polarized ground state | ft) 
of the Hamiltonian satisfying s + (x) \ ft) = for all x. 
In the spectral representation Eq. ([4]) takes the form 

S{k,w)=Y t S(w-Ef{k))\{f,k\a-(k)\ ft)| 2 , (5) 
/ 

where the sum is taken over the eigenstates \ f, k) of the 
Hamiltonian j3| carrying the momentum k. The ener- 
gies E f (k) are defined by H\f, k) = E f (k)\f, k). The fre- 
quency LU-(k) in Eq. JT]) is given by u>-(k) = minj Ef(k). 
Thus the calculation of u)-(k) reduces to the analysis of 
the energy spectrum of excitations. The calculation of 
A(fc) directly from the formula {5J is a far more difficult 
task. It requires the knowledge of the matrix element 
and their resummation procedure. For most integrablc 
models, including Yang-Gaudin, such calculation is be- 
yond the reach of the existing theory. A way to bypass 
this problem is to combine the BA with an effective field 
theory. This is the route we take in our calculations. 

We begin our analysis with a brief description of BA 
equations and a calculation of u>-(k). All the states \ f, k) 
in Eq. j5]) lie in the sector with the z projection of 
the total spin given by S z = N/2 — 1. In this sec- 
tor Bcthc's wave functions are characterized by a set 
of quasimomenta {Ai, . . . , Ajv, £,} which satisfy the BA 
equations [3] 



JV 



LXj + S(Xj -X k ) = ZttIj + 9{2Xj - 2f) 



(6) 



it=i 



Here 9(X) = 2 arctan(A/g) is the two-particle phase shift, 
and Ij = rij — (N + l)/2, where rij are a set of dis- 
tinct integers. The branch of 8(\) is chosen so that 
#(±oo) = ±7r. The total energy E and momentum P of 

a system are given by E — J2f=i an( l P = Yj4=i -\n 
respectively The quasimomcntum £ enters in E and P 
indirectly, through the solution of Eqs. ([6]). In the limit 
£ = oo Bethe's equations ((6]) are identical to Bethe's 
equations of the fully polarized system [2^|, S z = N/2, 
which is equivalent to the Lieb-Liniger model 2J, |25(. 



The distribution of Ij in the ground state of the model is 



3 



N+l 



3 = 1, 



, N. 



(7) 



Introducing the quasimomenta density p(Xj) = 
l/[L(Aj+i — Xj)] and taking the thermodynamic limit 
< po < oo as N, L — > oo one gets the integral equation 



(8) 



for the quasimomenta in the state (JTJ) and £ = oo. The 
kernel K(\v) = K(X - v) is K(X) = 86(X)/dX = 
2.9/(ff 2 + A 2 ). Note that p(X) should satisfy f* A dX p(X) = 
po- This formula together with Eq. ([5]) is used to get the 
value of the Fermi quasimomentum A as a function of 
the particle density po. The ground state energy in the 
thermodynamic limit is 



En =L 



dX X 2 p(X) 



(9) 



and the momentum of the ground state is zero. 

Consider now the state characterized by a finite value 
of £ and Ij given by their ground state values, Eq. ([7|. 
This state is an excitation above the vacuum, which we 
shall call a magnon. Introducing the so-called shift func- 
tion [111 by F(\j\£) = {Xj - Xj)/(X j+1 - Xj), where A^ 
are ground state quasimomenta, and Xj are those of the 
excited state, we get the following integral equation for 
F in the thermodynamic limit: 



tt + 9(2X - 2£) 



2tt 



(10) 
(11) 



The momentum of the excited state is 

k = [ dXp{X)[ir + 9{2X-2£)}, 

J-A 

and its energy above the ground state is 
1 f A 

u)_(k) = / dAe(A)X(2A-2£). (12) 

K J-A 

Here ui- is written as a function of the physical (observ- 
able) momentum k, which is related to the quasimomen- 
tum £ by the integral equation (jTTJl . The quasiencrgy 
e(A) is given by the solution of the integral equation 



•(A) 



2tt 



dv e(v) K(X, v) 



X 2 



P 



(13) 



-A 



satisfying a condition e(±A) = 0. The parameter p en- 
tering Eq. (fT3|) is the chemical potential, defined by 
p = (dE /dN) L , where Eq is found from Eq. ©. One 
can show that at small k the dispersion law (I12p is 
parabolic [27|, aj_(fc) = fc 2 /2m*, with the effective mass 
satisfying m" 1 = -(Trgpg)" 1 J A A dAe(A). 

Another way to excite the system is to create a particle- 
hole pair by moving one of the quantum numbers Ij in 
© outside of the ground state distribution ([7]) . Such ex- 
citations are analyzed in detail in [25| (see also Ref. [26|). 
In particular, at small momentum they are shown to be 
equivalent to sound waves propagating at velocity 



1 



de(X) 



2tt P {K) dX 



(14) 
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FIG. 1: Excitation spectrum in the model ([3]). Panels (a) and 
(b) show the dispersion curve, Eq. (|12|l . of the magnon 

(solid line) and the exact lower bound of the particle-hole 
continuum (dashed line) for two values of 7 = g/po- Panel (c) 
shows the 7 dependent exact lower bound E p of the particle- 
hole continuum at k = /cf = irpo- Panel (d) shows the 7 
dependent magnon energy E fd — LO— at k — kp. 



Any excitation in the ./V-particle sector with S z = 
N/2 — 1 consists of several particle-hole pairs and one 
magnon. The exact lower bound of the particle-hole con- 
tinuum [25| and the dispersion curve of the magnon are 
illustrated in Fig. [T] (a,b). For all values of coupling con- 
stant g the magnon branch lies below the particle-hole 
continuum. It is thus single magnon dispersion, Eq. (|12[) . 
which gives exact lower bound of the excitation spectrum. 

While w_(fc) is found using BA exclusively, in order 
to get the threshold exponent, A(fc), in Eq. (fTJ) for the 
function we need to combine the BA solution with a 
low-energy effective field theory. To do so, we employ the 
method proposed in [l7| . We introduce an auxiliary mi- 
croscopic theory with a local Hamiltonian H depending 
on k as an external parameter and having the following 
properties: (i) it conserves the total momentum, which 
will be denoted by q. (ii) its excitation spectrum at q = k 
is gaplcss. (iii) its structure factor S satisfies 



<%,w) 



S(q,to~(k) + oj) 



1. 



k, oj^O. 



(15) 



In integrablc models H can be constructed as a linear 
combination of a finite number of mutually commut- 
ing local integrals of motion. The eigenstates \f,q) of 
H are at the same time the eigenstates of H, there- 
fore S(q,uj) = Z f 6(u - E f (q))\(f,q\s.(q)\ 1\)\ 2 , where 
H\fiQ) = F/(g)|/, g). Like for H, the low energy spec- 
trum of H consists of sound waves and a magnon. The 
energy of the magnon is proportional to (k — q) 2 as q — > k. 
The condition (|15p requires that the velocities of the 
right- and left-moving sound waves be different and given 



v± = v s ± duj-(k)/dk, 



(16) 



where v s is given by Eq. (|14| . 

The dynamics of sound waves is governed by the Lut- 
tinger Hamiltonian 



J^Hr, H r = -^ dx: [d xVr (x)] 2 :, (17) 



where the operators <p r are chiral boson fields, 
[ip r (x), (p r i (x 1 )] = iTrr6 rr isgsi(x — x') related to the mi- 
croscopic particle density by 



p{x) = po + {2n)-WK[d x <p + (x) - d x tp-(x)] (18) 

and the symbol :: stands for the boson normal ordering. 
In order to describe the low-energy magnon excitation 
we introduce the spin density field s(x), related to the 
microscopic spin density by s z (x) = s z (x) + po/2 and 



s±(x) 



s±(x), where s± = s x ± is y are the lo- 



cal spin-ladder operators of Eq. ([4]). Within the effective 
theory the operators s± are smooth spin flip fields. Since 
a local spin flip may excite sound waves, an effective the- 
ory should contain a coupling between s and tp±. The 
minimal local coupling respecting the SU(2) symmetry 
of the microscopic theory and vanishing in the absence 
of magnon excitations 3l| is 



Hi 



E 

r=± 



v r /3 r 
~2^~ 



dx d x ip r (x)s z (x). (19) 



Other possible couplings involve higher gradient terms, 
which do not contribute to the critical exponents. The 
kinetic energy density of the spin field is represented 
by a higher gradient term d x s+(x)d x S-(x) that can 
also be neglected in the calculation of the critical ex- 
ponents [28j . The total Hamiltonian of the effective the- 
ory describing the dynamics near the threshold is thus 
given by H Q g — Hq + Hi. This Hamiltonian is diag 
onalized by a unitary transformation e 
5 = (27T)- 1 J Q L dx[/3 +( p+(x) - 
function (01 this gives 



s H e ffe~ iS with 
x)]s z (x). For the 



A(fc) 



1 



(20) 



What remains is to determine the coupling constants f}± 
in terms of the parameters of the microscopic theory. 
This is done by the comparison of the low-energy spec- 
trum of the microscopic Hamiltonian H , found from the 
BA solution, and the spectrum of the effective Hamilto- 
nian H c g. This procedure yields 



= 27rrF(rA|0, 



±1, 



(21) 



where F is defined by the solution of the integral equa- 
tion (|10|) . We solve this equation and find A(fc) nu- 
merically for different values of the coupling constant 



4 




FIG. 2: The function a(k) denned in Eq. (26) is plot- 
ted for different values of the dimensionless coupling con- 
stant 7. The values of the Luttinger parameter K are in- 
dicated for each curve and correspond in increasing order to 
7 = 00, 1.65, 0.56, 0.238 and 0.109 respectively. 



7 = g/po- For easier comparison with Eq. ^ we rep- 
resent our result in the form 



A(fc) 



K f k 



2 \k 



K 



-a(fc), (22) 



where kp = npo and K = kp/v s is the Luttinger pa- 
rameter calculated using Eq. (fl"4|) . The function a{k) for 
different values of the Luttinger parameter is shown in 
Fig. [2] One can see that at large values of 7 (or K close 
to one) the small k expansion derived in Ref. [B| is valid 
for all values of k, in agreement with Ref. Q- It is inter- 
esting to note that the leading correction to this result is 
of the second order, (K — l) 2 ~ j~ 2 . At arbitrary inter- 
action strength a(k) ~ fc 4 as k — > (^.therefore the small 
k expansion of A(fc) found in Ref. [5j remains valid for 
all values of 7, confirming the general result of Ref. [1]. 
Note that a(k) also vanishes at k = 2kp. 

The problem considered in the present work is directly 
related to the X-ray edge problem in the theory of the 
mobile impurity. In this context, the model ([3]) was in- 
vestigated in Ref. [22]]. The approach of Ref. [22| ex- 
ploits a transformation to the co-moving reference frame 
and combines BA with an effective field theory similar 
to ours. The method of Ref. 22j has recently been suc- 
cessfully applied to the Heisenberg model and later was 



shown 19 1 to produce results equivalent to the method 
of Ref. 17| used here. A direct comparison of the present 
work with Ref. [22| is however not possible, because the 
latter used an incorrect BA solution of the model j3]) . 
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